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Abstract 

We construct a microscopic model of electron interactions which gives rise to both 
superconductivity and antiferromagnetism, and which admits an approximate SO{b) sym- 
metry that relates these two phases. The symmetry can be exact, or it may exist only in 
the long-wavelength limit, depending on the detailed form of the interactions. We compute 
the macroscopic Landau-Ginzburg free energy for this model as a function of temperature 
and doping, by explicitly integrating out the fermions. We find that the resulting phase 
diagram can resemble that observed for the cuprates, with the antiferromagnetism real- 
ized as a spin density wave, whose wavelength might be incommensurate with the lattice 
spacing away from half filling. 



1. Introduction and Summary 



Zhang has recently proposed [1] that the superconducting and antiferromagnetic 
phases of the high-Tc cuprates might be related to one another by an approximate 5*0(5) 
symmetry of the electronic Hubbard Hamiltonian. This proposal incorporates in a funda- 
mental way the connection between antiferromagnetism (AF) and superconductivity (SC) 
which is observed in these systems. 

Recently, however, the foundations of this picture have come under attack, with crit- 
icism directed against the existence of an approximate SO (5) symmetry of the Hubbard 
Hamiltonian [2] (for a reply, see [3]), as well as against the very possibility, in princi- 
ple, of relating the antiferromagnetic and superconducting phases by the rotation of a 
finite-dimensional order parameter [4]. This motivates studying whether and under what 
conditions an approximate SO{5) symmetry can emerge from a microscopic picture of 
electron dynamics. 

In this paper we investigate a model of electron dynamics which can naturally incor- 
porate an approximate SO (5) symmetry.^ Our model consists of degenerate electrons with 
two kinds of attractive interactions, described by a four-Fermi Hamiltonian of the form 
H = Ho + Hint, where i^int = H^j, + Hgc, and 

^o = Y^ {Sp - IJ) i^li^p + ^ (a% + a° n* • , (1) 

p q 

Haf = -^^f{P^ Q) (i^l+Q+q^Ipp) ■ flQ+q + h.C, (2) 

pq 

Hsc ^^'^g{p,q) (^V-p-gf^2V'p)0Q + h.c., (3) 



2V- 

pq 



Here ipp — { ) is the electron field, while 0q and fiq are Hubbard-Stratonovich auxiliary 



scalar fields which, when integrated out, produce the two four-Fermi interactions. Sp 
denotes the electron single-particle dispersion relation, fi is the chemical potential that 
measures deviations from half-filling, and a — (ui, i — 1,2, 3} are the Pauli matrices. Q is 



^ A broader class of similar models has recently been proposed in ref. [5] . 



a fixed nesting vector which will be defined below, and a° are positive constants which 
could, if one wished, be absorbed into the definitions of / and g hj a redefinition of (f) and 
n; however it is convenient to keep them for the purposes of renormalization in §4. 

Our goal in this work is twofold: first, to find the conditions under which this model 
exhibits an approximate SO{5) symmetry, and to see how the symmetry-breaking effects 
manifest themselves at temperatures much less than the Fermi energy. We find the sym- 
metry to be possible, and it can arise "accidentally" at long wavelengths even if it is not 
important for the microscopic dynamics. We also find that symmetry-breaking effects are 
marginal, for weak couplings, in the sense that they increase only logarithmically as one 
scales into the infrared. Second, we wish to compute the effective Landau-Ginzburg (LG) 
potential, so as to determine how the phase diagram depends on the microscopic electronic 
couplings. The result is that the AF phase can have the phenomenologically desired shape 
in the temperature-doping plane, preferring low temperatures and zero doping. Moreover 
the SC phase, while also preferring small doping, is not suppressed as quickly by increased 
doping as is the AF phase. The result is that SC becomes energetically favorable to AF 
at some lower critical doping, and persists until some upper critical doping. These results 
are consistent with but more predictive than the general effective-field-theory description 
of [1] and [6]. For example they predict the absence of a mixed AF/SC phase away from 
zero doping when the fundamental Hamiltonian is /SO (5) -invariant. They also agree with 
the observed phase diagram of the high- Tc systems, although we do not attempt here a 
detailed discussion of the applicability of these models to the cuprates. 

Our analysis of the model defined by eqs. (l)-(3) is motivated by the renormalization- 
group (RG) approach to understanding superconductivity and antiferromagnetism as BCS 
and spin-density- wave instabilities within Fermi liquid theory [7]. In the RG language, 
a Fermi liquid is understood as a regime of scales for which the dominant quasiparticle 
degrees of freedom are weakly-coupled, degenerate fermions carrying the same quantum 
numbers as the underlying electrons. Under such circumstances almost all of the quasipar- 
ticle self-interactions are irrelevant in the RG sense: they become less and less important 
as one integrates out high energy modes to obtain an effective Hamiltonian valid near the 
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Fermi surface. 

The only exceptions to the rule that interactions are irrelevant in the infrared are 
certain four-Fermi terms, which can be marginally relevant for special kinds of electron 
kinematics. An important example is any pairwise attraction between electrons (or holes) 
having opposite momenta. This kind of interaction grows logarithmically in the infrared, 
eventually becoming strong enough to trigger the BCS instability at sufficiently low ener- 
gies. A second exception exists when the quasiparticle Fermi surface (FS) is nested.^ In 
this case the attraction between electrons whose momenta sum to the nesting vector, Q, 
is also marginally relevant, potentially triggering an instability towards the formation of 
a condensate which is modulated in space with wavevector Q. The result is a charge- or 
spin-density wave depending on the electric charge and spin of the attractive channel. 

We have both types of instability in mind when using the Hamiltonian in eqs. (l)-(3). 
It is assumed that the quasiparticle energies Sp are time-reversal invariant, £_p = £p, and 
have a FS with a nesting vector Q, defined by the property £p+Q = —Sp for p near the FS. 
(Our convention is to define the FS at half filling to be the zero of energy.) An example 
of a dispersion relation with these properties in two dimensions is that of the extended 
Hubbard model. 



general analysis that follows in §4 does not depend on using this explicit form. Notice that 
because antiferromagnetism arises here as a spin-density wave, its periodicity need not be 
precisely commensurate with the lattice spacing, particularly away from half filling. Such 
incommensurate antiferromagnetism is particularly interesting in view of the splitting of 
the AF peaks which is seen in recent neutron-scattering data [8] . 

Since eqs. (l)-(3) involve an attractive interaction in the spin-triplet channel, the nest- 
ing produces a spin-density- wave instability. The Hubbard-Stratonovich formulation of the 
Hamiltonian, eqs. (l)-(3), is particularly convenient for following these instabilities of the 

The FS is nested if its opposite edges are related to one another by a fixed tranlation, Q, in momentum 




(4) 



where a is the lattice spacing and the nesting vector is Q = (± ^, ± ^) . However the 



space. 
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system since they may be studied by examining the minima of the potential for the modes 
00 and ng, which is induced when the electrons are integrated out. Depending on which 

of these modes condenses most strongly, the resulting ordered phases are superconductors 
or antiferromagnets, possibly both. 

Our presentation is organized as follows. In the next section, §2, the conditions for 
the Hamiltonian, eqs. (l)-(3), to exhibit SO{5) symmetry are derived. We find that the 
symmetry is present in the long-wavelength limit, g — > 0, if the coupling functions satisfy 
the relation p = . In §3 we integrate out the electrons to obtain the long- wavelength 
effective potential for the order parameters, 0o and uq. The expressions obtained in 
§3 are evaluated in §4, in the limit of degenerate electrons, i.e. with all couplings and 
dispersion relations linearized around the FS. We present explicit expressions in this limit 
for the coefficients of the Landau- Ginzburg theory, as functions of the fundamental fermion 
couplings, temperature and chemical potential. In some parts of the phase diagram the 
usual expansion of the free energy to quartic order in the fields is unbounded from below, 
forcing us to use the full potential. From these results it is straightforward to determine 
the main features of the phase diagram. We find the dependence on T and \x can have 
the desired form if the AF interaction in eq. (2) is chosen to be somewhat (~ 10 — 20%) 
stronger than the SC interaction of eq. (3). In §5 we examine the validity of the degenerate- 
electron approximation by recomputing the quadratic and quartic terms of the effective 
potential without linearizing about the FS. We do so using specific choices of e, / and g 
that have been previously suggested [1,5,9] in connection with d-wave superconductivity, 
to show that the approximations used in §4 are indeed reliable. §6 concludes with a brief 
recapitulation of our results. 
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2. The SO (5) Symmetry 



In this section we identify the conditions under which the Hamiltonian, eqs. (l)-(3), 
admits an SO{5) symmetry along the hnes Zhang has proposed. 

2.1) The Algebra 

SO{5) symmetry contains as a subgroup the S0{3) x S0{2) group of spin rotations 
and electromagnetic gauge transformations, with respect to which H is obviously invariant. 

The electron part of the corresponding conserved charges (electric charge and spin) are 
given by: 

p p 

The nontrivial requirement for SO {5) invariance is the existence of the electrically- 
charged, spin-triplet off-diagonal operators, H, whose real and imaginary parts fill out the 
generators of SO (5) in addition to Q and S. Motivated by Zhang's [1] choice, we take its 
action on electrons to be given by 

n = ^ /ip ('i/'-p+Q(^2 O^-iAp), (6) 

p 

where hp is a function to be determined. As is easily verified, this has the desired commu- 



tation relations with Q and S, and satisfies 



na,n6 



= 0. In addition, if the function hr 



shares the same properties as does Sp: i.e. hp+q = —hp, and h-p — hp, then: 

W^cV'p)]- (7) 



na,n,* 



Eq. (7) shows that {Q, S, 11,11*} satisfy the SO{5) commutation relations provided we 
choose \hp\ = c, where c is an appropriately chosen constant. This is only consistent with 
the property hp+Q = —hp if hp = ±c for p inside the FS, and hp = outside. 
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2.2) Conditions for Symmetry 



We next determine the conditions under which the Hamiltonian, eqs. (l)-(3), com- 
mutes with the n. A nonvanishing chemical potential explicitly breaks 5*0(5) down to 
SO{3) X S0{2), so we set = here. Furthermore, we ignore for now the scalar part of 
Hq, which is invariant provided and n transform together in the fundamental represen- 
tation of SO{5). The electron part of the free Hamiltonian then satisfies: 



-2 Sp hp ilpQ_p (J2(Ja ^p)- 



(8) 



This vanishes by virtue of the Fermi statistics of V', provided that Ep shares the previously 
advertised symmetry properties of hp-, i.e. if Sp+g = —Sp and S-p = Sp. 

For notational convenience, we define the operators appearing in i^int by Vb{q) and 
W{q), where Vb{q) = T.pf{p,q) {ipl+q+Q CTb i^p) , and W{q) = Ep^(P, (V'-p+gf^2 V'p)- 
Then 



^a,W*iq) 
^a,Vk{q)\ 



= '^Y[9*{p + Q - Q,q) + 9*{Q - P,q)\ hp (V'p-g+g^^aV'p), 

p 

= 2 ^ fip, q) hp+q (V'Ip-q 0-2 CTaO-b V'p) 

p 

= [j'^P' ^) + ^^~P ~ ^' ^) ^-p] i'^-p-q ^2 V'p) 

p 

+ i^abc Y l^^P' ^) ~ ^^~P ~ ^' ^) ^-P i'^-p-q ^2 CTc V'p)- 



(9) 



The operators Va and W rotate into one another as a five-dimensional vector of 50(5) 
provided the coupling functions satisfy the conditions f{p,q) hp+q = f{—p — q,q)h-p oc 
9*{P, q) and [g*{p + q - Q,q) + g*{Q - p, q)]hp ex /(p, q). 

Of particular interest for many purposes, particularly the phase diagram, is the 
long- wavelength limit, g ^ 0. In this limit we have Ha, V6(0) = —A5abW*{Q) and 
Ila, (0) = -4 Ki(0), provided the functions fp = limg^o fip, q) and Qp = limg^o 9{p, q) 
satisfy the properties fp = f-p, Qp = g-p = -Qp+q and fphp = Qp and Qphp = fp. Given 



our earlier choice that hp = ±c, with opposite signs inside and outside the FS, we see 
that W{0) and Vb{0) fill out the fundamental representation of SO (5) provided that the 
coupling functions satisfy fp — sign{e p)gp. 

Finally, so long as W and Vb transform as a 5 of S0{5), the entire Hamiltonian, 
eqs. (l)-(3), is SO{5) invariant, provided that S0{5) is taken also to rotate the five scalars, 
Re{(f)o) I a/2, Im((/)o) / \/2 and ng, into one another as a 5. Notice that although the demands 
made on / and g for exact S'0(5) symmetry are quite strong, the symmetry can emerge 
for small q under much weaker conditions, and this is sufficient for the invar iance of the 
phase diagram. 

3. Integrating out the Electrons 

We now turn to the derivation of the effective theory which governs the long- wavelength 
scalar modes, 0o and uq. The resulting effective potential will be subsequently used to 
identify the phase diagram of the system, and to see to what extent it exhibits SO{Ji) 
invariance. Since the effective potential must be symmetric under IJ{\) transformations of 
00 and under 5'0(3) rotations of ng, for notational simplicity we shall take 0o = real 
and {nQ)i = nSi^^. 

3.1) Preliminary Issues 

The thermodynamic variables we wish to follow in the LG potential are the values of 
the fields and n themselves, as well as the temperature, T, and the electron doping away 
from half-filling, x. We therefore ignore all scalar modes apart from (j) and n, and compute 
the effective potential for these variables. In order to follow the doping we also introduce 
a chemical potential, /j,, which measures the change in density away from the free-particle 
FS (which is defined here by Ep = 0). 

We now integrate out the electrons, using the Hamiltonian, eqs. (l)-(3), supplemented 
by the chemical potential. To incorporate the Fermi sea we exchange the electron field, ipp 
for p inside the FS, for a hole field, Xp = (^2'^-pi in the usual way, permitting the use of 
vacuum propagators for the integration over the fermion fields. 
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For the purposes of performing the functional integral over ipp and Xpi it is convenient 
to group 'i/'p, V'p and their counterparts having momentum —p and p — Q into the following 
16-component vector: 







( ^P \ 


i^p-Q 




^2 XQ-p 


%-Q 




(72 XQ-p 






V'-p 


V'-p 




r-p 


i^Q-p 




^2 Xl-Q 




\ (72 Xp-Q ' 



(10) 



When '^p is used as the field, all sums over p must be restricted to the original Brillouin 
zone, modulo the transformations p —p and p ^ p + Q. 

With this choice the quadratic part of the Hamiltonian may be written | ^j^A^p, 



with: 



A = 



A B 
B A 



and A and B denoting the following 8x8 matrices: 



(11) 



A 



fSp-jJ, fpU (73 \ 

-Ep + 11 fpU (73 



fpn (73 

V 





fpn (73 



-£p - H 









(12) 



B 



(13) 



and 

/ gp(t)(T2 Q \ 

9p<p(72 

gp(j)a2 

V gp(j)a2 / 

In these last expressions, each entry of the matrices A and B are 2x2 matrices in spin 
space. 



3.2) The Functional Determinant 



Standard methods may now be used to compute the functional integral over the 
fermion fields. Using the matrix form for the Hamiltonian, eqs. (11)-(13), the result 
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for the one-foop, finite-temperature contribution to the free energy density, Fi, in d spatial 
dimensions is: 

-PFi (0, n) = ^ In det [-0^ + 

4/|f.EK-l+^i)+>-H+^^-)]' (14) 



3 = -oo 

In [cosh {(3 A+/2) cosh {(3 A_/2)] , 



(27r) 



d 



where, as usual, (3 = 1/kT and the sum is over the Matsubara frequencies, Uj = (2j + l)7r/?. 
We have switched to continuum momenta, and restored the integration region to the full 
Brillouin zone (so that it is no longer modded out by p — > —p and p P + Q)- The 
eigenvalues of the Hamiltonian A are given by 



A± 



2 i 

{^/4+{M^±f^) +{9p4>fy- (15) 



Even before undertaking a detailed evaluation of the free energy (14), we can prove 
an interesting result: in the SO(5) limit where fp/a^ = Qp/a'^, and for nonzero doping, 
/U 7^ 0, the ground state must have ^ = or n = (possibly both) and hence is never 
a mixed phase of AF and SC. This follows from the form of the energy eigenvalues (15) 
in the conditions for minimizing the free energy. Adding Fi to the tree-level free energy 
Fq = a^cfy^ + a°?7,^, these conditions are 

- =a%X-- —^y tanh(/?A±/2)— ^ (16) 

where X stands for either of the fields or n. The derivatives of the eigenvalues have the 
form 



aA±^^_ d\± ^fpn 1^ fi 



(17) 

When fp = Qp and a^^ — a^ — aP, the minimization conditions can be written in the form 
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where it is easily shown that G and jiH are positive definite functions. Therefore in this 
limit there are never simultaneous solutions where both and n are nonvanishing, except 
when jjL vanishes. 

When evaluating these expressions in the next section, we focus on the degenerate 
limit, for which kT is much less than the Fermi energy at zero doping. For example, using 
the explicit dispersion relation given by eq. (4), this corresponds to neglecting powers of 
kT /t. In this limit only states near the Fermi surface are important in the integrals, and it 
is useful to decompose the momentum p into a component k on the FS and a component 
I perpendicular to the FS. We may then linearize all quantities in restrict the coupling 
functions g and / to momenta lying in the FS, and obtain for the dispersion relation 
£ ~ Vp(., with the Fermi velocity, Vp^k), generally being a function of momentum, k. 

Eqs. (14) and (15) are this section's main results. We now extract their implications 
for the LG potential,^ first by examining them in detail in the vicinity of zero fields, 
</) = n = 0, and then by computing the potential for arbitrarily values of the fields, in 
certain directions in the c()-n plane. 



4. The Landau-Ginzburg Theory 



Eqs. (14) and (15) define the effective potential for the SC order parameter, 0, and 
the AF order parameter, n. Since this potential is analytic^ at = n = 0, an expansion to 
quartic order in the fields gives a good indication of when a condensate will form for either 
(/) or n. However we will find that the quartic terms do not always stabilize the potential, 
and it is necessary to know how the full potential looks away from cj) — n — In this case 
we will use the full free energy in order to determine the phase diagram of the system. 



See refs. [9] for related calculations of the LG free energy for high-Te systems. 

The same is not true of the zero-temperature limit of the potential, which has the nonanalytic behavior 
In (jp in the vicinity of <^=0. 
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4.1) Analytic Expressions in the Quartic Approximation 



We start by approximating the free energy with a quartic potential, whose form we 
take to be 

(3[F - F(0)] - (a^ (I? + a^n^) + U^^^^^ + 2 h^n4^''n' + hnnU^) + • • • (19) 

where the eUipsis denotes terms involving more than four powers of (p or n. Using Eqs. (14) 
and (15) and defining 



Nn 



_ 1 f d'^'krg- tanh(/3e/2) 
'^"^y (27r)V(/c)' e ' ^ ^ 



we find that the quadratic coefficients are given by the following integrals 

fA 



= al{A) - No,2 f deTp{0, (21) 
an = a° (A) - N2,o t de T^(0-, (22) 

J-A £ 



where and are the "bare" values of these quantities that appear in the original 
Hamiltonian, eqs. (l)-(3). The integral over e diverges logarithmically in the ultraviolet, 
forcing the introduction of the cutoff scale, A. Physically this scale represents the energies 
for which the linearization of the quasiparticle spectrum no longer applies, or where the 
degrees of freedom of the Hamiltonian, eqs. (l)-(3), are no longer appropriate, whichever 
is lower. If these approximations were exact, then A would be the actual maximum energy 
available on the lattice. In any case, we can consider A to be of order the Fermi energy. 
The divergent dependence on A as A — > oo can be absorbed into a renormalization of 
and a°. 

The quartic coefficients are given similarly by the following expressions, which since 
they are convergent allow us to take the cutoff A to infinity if we wish: 

/A J" 
de^^, (23) 
-A ? 
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/A J"' rd\ 
de (24) 

.„„^-../>(?M-f|(^)). 

The approximation A — > oo is only valid if T and fj, are much smaller than the Fermi 
energy. 

Because we work in the degenerate limit, near the FS, and at zero field, the shape 
of the coefficient functions as k varies over the FS enters only as an overall normalization 
factor. As expected from our previous discussion, the /SO (5) symmetry is manisfestly 
present in the limit = and /^/a° = g^/a^. 

4.2) Instability Towards Condensation 

Before turning to the numerical evaluation of these integrals, it is instructive to first 
explore the physical implications of the renormalization itself. Because the sign of the 
loop contributions to a<^ and is negative, these quantities decrease as the temperature 
is decreased. In the limit of large /3A, it is easy to see that the one-loop contributions 
Agx — ttx — have the limiting behavior 

Aa^ ~ -2Aro,2ln(/?A); Aa^ ~ -2N2,oHP^)- (26) 

Even if a^(T ~ A) and a^iT ~ A) are initially chosen to be positive, so that F is minimized 
hy (p — n — 0, eventually their decrease with decreasing T can drive them negative, causing 
an instability towards the development of nonzero condenstates for or n. 

Eq. (26) gives the evolution of the 50 (5) -breaking combination of the electron cou- 
plings: 

pQ-piH - an) = 2iN,,o - No,,) = 2 j ^-^^ (f - g% (27) 

This shows that any 50(5) breaking present in the original Hamiltonian due to p ^ 
will be enhanced in the effective GL Hamiltonian by a potentially large logarithm of the 
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form ln(/?A), where A is of the order of the Fermi energy. For example, suppose the AF 
pairing is taken to be s-wave, / = /o, whereas the SC pairing is d-wave, corresponding 
to g = go{cos{akj:) — cos{aky)). Here /o and go are arbitrary constants. In general A''2,o 
and A^o,2 differ from each other, and although go/ fo can be tuned to make the difference 
vanish, there seems to be no symmetry principle that would require it to take this special 
value. Moreover, even if such a choice were made, it would still not insure that the quartic 
couplings would have SO{5) symmetry, i.e., A^4,o = -/Vo,4 = N2,2- In §5 we will say more 
about the extent of S'0(5)-breaking that arises from this choice of coupling functions. 

4.3) Numerical Results near = n = 

Figures (1) and (3) present plots of the coefficients of the potential as functions of 
the chemical potential, factoring out the integrals over the FS. In this section we describe 
these plots in more detail, and use them to deduce the general shape of the phase diagram 
in the temperature-doping plane. 

Figure (1) plots the one- loop contribution to the quadratic coefficients in the dimen- 
sionless form (a^ — a°)/A^o,2 and (a<^ — a°)/A^2,0) where A^n,m are the FS integrals defined 
in (20). At zero doping these functions are negative and are increasing with temperature 
like ln(T/A), so it is always possible to choose the bare coefficients or in such a way 
that a„ or is positive for high temperatures and changes sign at some Tg, triggering 
the instability toward condensation. As the doping is increased, increases much more 
rapidly than a,^. It follows that if the AF phase is energetically preferable to the SC phase 
at // = 0, the system will be AF at zero doping, and then make the transition to SC at 
some small doping, provided that is still negative at this value of 

With applications to the high-Tg cuprates in mind, we imagine A ~ 0.1 eV ~ 10^ K 
and so we adjust the bare coupling a° to ensure that the AF transition (Neel) temperature, 
Tjv, at fi — is 1/10 of the cutoff scale, /9jvA = 10. Prom fig. 1 it can be seen that at a 
temperature T = T^/10, for example, we have /3A = 100 and so the AF phase is quenched 
at a critical doping of Xc ~ A*c/-^ = 10%. This is the doping at which a„(//c )/3A=ioo 
becomes equal to an(0)/3A=io- 
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Figure 1 

The coefficients (a,^— o^)/iVo,2 (dashed) and (a„— a° )/iV2,o (solid) as a function of for /3A=10, 

20, 50, 100 and 200. 

The above example ignores the SC order parameter, which would be valid if was 
taken to be much larger than a^. If on the other hand SO{b) symmetry was realized, so that 
aP^ = a° and A^o,2 = -^2,0) then at /x = the ground state would be an arbitrary admixture 
of AF and SC phases since both have the same energy. In this case, the system will always 
make a transition to the pure SC phase at the critical doping Xc- Since experiments 
indicate that at = the cuprates are pure AF, the SC phase must be suppressed 
by having a^/A^o,2 > an/^2,05 so that SC is energetically disfavored until the critical 
doping is reached. In this situation, the phase diagram must have a form similar to figure 
2. Of course, this prediction does not take account of quantum fluctuations within the 
effective theory itself. Such higher-loop effects are believed to be especially important 
near the bicritical point, where the SC phase first emerges [1]. They may have the effect of 
driving this point downward, leaving a gap between the SC and AF phases at temperatures 
somewhat above the existence of the bicritical point. 

To determine exactly how much 5*0(5) breaking is needed in the quadratic terms in 
order to suppress SC at zero doping, one must compare the local minima of the potential 
along the = and n = axes and insure that the minimum with n 7^ is the deeper 
of the two. This issue can usually be settled by looking at the higher order terms in the 

15 



Figure 2 

Schematic representation of the phase diagram for the present Hamiltonian, with SO (5) symmetry 
broken to favor the AF condensate. Quantum effects not considered here may be able to drive the bicritical 
point downward. 

potential. To tliis end we liave computed tlie quartic coeflacients in tlie dimensionless forms 
b^^T^/No^4, bnnT^/N^,o and bn4>T^/N2,2, and plotted them in fig. 3. 

The fact that 6,^0 is only weakly dependent on the doping can be see by observing 
that when A — * cxd, // can be completely removed by shifting the integration variable from 
£ to ^. The low temperature form of the quartic couplings can be further understood 
analytically through the identities 



lim 



lim 



7C(3) 



27r2 



(28) 
(29) 



Eq. (29) is valid so long as it is combined with a function that is nonsingular at £ = 0, in 
which case its own 1/e singularity can be integrated over using the principal value. Both 
of these representations of the delta function can break down if multiplied by singular 
functions of ^, which is precisely what happens at = 0, but for ji^l/ (3 they can be used 
in conjunction with (23)— (25) to deduce that hnn = bncj) = and 60^/Ar4,o - 7C(3)/27r2, in 
the low temperature limit. 

At zero doping, all three coefficients have essentially the same value, and from the 
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preceding paragraph we understand why b,^,^ remains nearly constant up to very large 
doping, whereas bnn and drop quite sharply. In fact, they reach negative values for 
intermediate dopings. For these dopings, the truncation to fourth order in the fields gives 
a potential which is unbounded from below. Thus the expansion of the potential cannot 
be used to determine the minimum in the n direction, and our goal of finding the absolute 
minimum cannot be attained by expanding only to quartic order. We must use the full 
functional form of the free energy for this purpose. 
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The coefficients 6^^T^/Aro,4, 6„„T^/iV4_o and 6„^T^/Ar2,2 as a function of fx/A, for A/T=10, 20, 50, 
100 and 200. The curves approach an asymptotic limit as A/T— *-oo. 



4.4) Numerical Results for Arbitrary and n 

In computing the full expression for Fi in eq. (14), we can no longer factorize the 
momentum integral into a product of the form J d'^^k G{k) J dsH{s) because the integrand 
is not a simple polynomial in fp and Qp. However we can still treat it similarly to our 
previous computations if we consider the special case where p and g"^ are constant. We 
do not expect this assumption to radically change the features of the more general case, 
since the /c-dependence of / and g has played no essential role in the preceding discussion. 
Constant g"^ still describes d-wave SC and is consistent with 50(5) symmetry, so long as 
g changes sign across the nesting surface, e.g. gp = sign(ep). 
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With this assumption, the free energy reduces to a single integral over the energy. 
We have evaluated it numerically for the relevant ranges of doping, temperature, and 
order parameters. The results are most conveniently expressed in terms of the following 
variables, for which the cutoff. A, is scaled out: 

4) = gc^jk; h = fn/A; x = ///A; 
a^ = al/{mo,og^y, an = a^J{mo,ofy, p = pA/2 

i = e/A; \±^{{^/i^T¥±x)^ + 4>Y^^; 

^F = ll\eJ cosh(/3A,)cosh(/3A-) \ 
(3 Jo \ cosh(/3(£ + x)) cosh(/3(£ — a;))y 

In the last line we subtracted a constant to make AF vanish at = n = 0. Then the full 
effective potential can be written in terms of 

F = F/(4iVo,oA2) = d^4>^ + dnh^ + AF{4>, h, /3, x). (34) 

In figure 4 we show F for six values of the doping parameter, x = 0, 0.2, 0.4, 0.6, 0.8 
and 1 (horizontal rows) , and along three directions in the (f)-n plane (vertical columns) . F 
is plotted as a function of (jP, {(jP + n^)/2, and n^, in the three respective columns. Each 
graph is shown for the range of temperatures given by /3A = 4, 10, 20, and 200, with F 
always increasing with temperature. We have chosen = 2 and cifj, = 1.6, hence 5*0(5)- 
breaking at the 20% level, in order to illustrate the phenomenologically desired dominance 
of the AF phase at low doping. At a; = 0, the lowest minimum clearly occurs for n 0. 
At a; = 0.2, the minima become degenerate and the transition to the SC phase begins. 
Further increasing the doping quickly removes the AF instability, and gradually quenches 
the SC phase as well, depending on the temperature. This is in complete agreement with 
the picture obtained from our analysis of the quadratic coefficients and confirms the phase 
diagram shown in fig. 2. 
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Figure 4 

The rescaled effective potential F for six values of the doping parameter, x (horizontal rows) and 
along three directions in the </)-n plane (vertical columns). F is plotted as a function of 4>^ , {4>^+n'^)/2, 
and n^, in the three respective columns. Each graph is shown for the range of temperatures given by 
/3A=4, 10, 20, and 200, with F always increasing with temperature. 
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4.5) Terms with Derivatives 



So far we have presented the one-loop free energy only for spatially constant order 
parameters. If c/){x) and n{x) are not constant, there will be additional terms depending 
on their spatial derivatives. Of these, the most interesting are the (V0)^ and (Vn)^ terms, 
since they are needed to compute quantum corrections coming from within the effective 
long- wavelength theory (two-loop and higher corrections). 

By calculating the relevant Feynman diagrams, it is straightforward to show that the 
quadratic terms in Fi for (pq and Uq with an arbitrary momentum q can be obtained from 
our previous formulas (21)-(22) by making the replacements 



where = Sp — and p is the loop momentum to be integrated over. Expanding in 
the external momentum q gives the expansion in derivatives of the fields. From this 
prescription it is clear that 5*0(5) symmetry will be realized in the derivative terms so 
long as f^{p, q) = g^{p, q) and = 0, since the symmetry is then manifest in the energy 
eigenvalues \± that appear in (15). 

The resulting integral over p depends on the exact form of the dispersion relation Sp 
and the coupling functions /p, Qp, so we will defer further evaluation until the next section 
where we consider a specific choice of these functions. 



Our analysis so far has relied upon the approximation that the most important con- 
tributions to the effective potential are coming from electronic excitations near the Fermi 
surface. This approximation presupposes the Fermi energy to be much greater than the 
temperature and the chemical potential. On the other hand this hierarchy of scales may 
only be a factor of 10 in the actual cuprate systems we are interested in. It is there- 
fore interesting to corroborate our conclusions without assuming that Fi is dominated by 




(35) 



(36) 



5. Exact Results in a Specific Model 
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the contributions from near the FS. In this section we will exactly evaluate the first few 
terms in the small-field expansion of F using a semi-realistic lattice dispersion relation and 
pairing functions, which have been previously suggested in the literature [1,5,9]. It will 
be shown that the results are in good qualitative agreement with those based on the RG 
approach. 

5.1) Method of Solution 

For our specific model, we take the lattice dispersion relation (4) of the extended 
Hubbard model and choose 

f(p, q) = 1, gip - 9/2, q) = (cosiapa;) - cos{apy))/2 (37) 

as would be appropriate for s-wave AF pairing and d-wave SC pairing. The factor of 2 in 
the definition of g is chosen such as to maximize the SO{5) symmetry near the undoped 
FS, as is shown below, and the q dependence is chosen to simplify the derivative terms 
to be computed at the end of this section. We further assume that the system is two 
dimensional, d = 2. 

The choice (37) is motivated by two considerations: (1) we wish to examine to what 
extent SO (5) symmetry survives when f"^ ^ g'^^ since it has been suggested that SO (5) can 
be a good approximate symmetry in the long- wavelength limit, even if it is not exact in 
the fundamental Hamiltonian; and (2) one would like to see how much the general features 
of the phase diagram depend on having exact SO (5) symmetry in the Hamiltonian. We 
will show that making g^ breaks the symmetry significantly due to fluctuations far 
from the Fermi surface, even though those close to the FS give contributions to the free 
energy which are approximately 5'0(5)-symmetric. Nevertheless, the phase diagram is not 
affected in a qualitative way, and should still have the shape shown in figure 2. 

The effective potential at g = is still given by the general expression (14), but now 
the momentum integrals range over the Brillouin zone, —iz/a < Px.,Py < / a, and there is 
no need to introduce any cutoff because the energy is explicitly bounded by —At < e < At. 
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As before, we expand F to quadratic and quartic order in the fields. Defining ^p = Sp — fx, 
this gives 



(2 



(38) 

Next one would like to reexpress the momentum integrals as an integral over the energy 
times an integral over the Fermi surface. This can be accomplished using the identity 



where u = cos{apj^) and v = cos{apy). The integrals over u and v can be done exactly, 
leaving an energy integral times a density of states function, 

/ d^pgi; = 4t dsJ\fm{e/At). (40) 

The density of states Mm is given in terms of complete elliptic functions K and E-^ ^ 
defining e = e/4i and o: = (1 + |£|)/(1 — |£|), the first few are 

N,{e)=-^K{l/a) (41) 
1 + \e\ 

M2{e) = (1 + \e\){K{\/a) - E{l/a)) (42) 

A4(£) = \oi{l - \e\f{{2a' + l)K{l/a) - 2{a^ + l)E{l/a)) (43) 

These are shown in fig. 5. Because of our choice of normalization for gp in (37), they 
all have the same asymptotic behavior near the undoped FS, £ = 0: jV"^ ~ ln(l/|2;|), 
independent of m. However they quickly diverge from each other for energies above or 
below the FS. This is interesting from the point of view of the 5*0(5) symmetry, because 
it shows that even though ^r^ 7^ 1, as is necessary to have exact SO (5), close to the FS 
the symmetry is nevertheless approximately realized. However the logarithmic singularity 



^ as they are defined in all standard mathematical references except for Abramovitz and Stegun and Math- 
ematica, which use a nonstandard convention. See for example [10] 
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The density of states A/q, M2 and M/^ as a function of £ = sj^t. 



near £ = in Mm is too weak to compensate for the S'0(5)-breaking differences between 
A/o and Ml away from the FS, as we will presently show. 

b.2) Expansion near (f) = n = 

Now the expressions for the coefficient functions in the small-field expansion of Fi can 
be written in a way which looks very similar to our previous results (21)-(25). Let us first 
define 



where a is the lattice spacing, and m is the quasiparticle mass in the vicinity of the FS, 
where the dispersion relation can be approximated as e = /2m. We then obtain 



4 2m 



(44) 



K = 



2' 




(45) 



(46) 



(47) 
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b(j}n — 



K 



K 



At rp, j'^N 

-u ^ 

-u \ ^ £ ds 



MO 



(48) 
(49) 



The dependence of these coefficients on doping and temperature is shown in Figures 6 
and 7, which can be compared to Figures 1 and 3 from our model-independent analysis. In 
Figure 6, which shows the 1-loop contributions to the quadratic coefficients, one sees that 
our previous result for an is in good agreement with the exact computation in the present 
model. This is not suprising because fp is treated as a constant in both cases, and the 
only difference between the two is the lattice dispersion relation in the present computation 
versus the approximation e ~ Vpl in §4. However the difference between the two approaches 
is quite apparent in a^, where Qp is relevant: increases with doping faster in the exact 
computation than in the approximation of factoring out the FS integration. Nevertheless 
the important qualitative feature needed to obtain the phase diagram of Figure 2 persists: 
ttc^ increases more slowly as a function of doping than does an- 

Similar remarks apply in comparing the quartic coefficients in the two treatments. 
The difference in the dispersion relation causes a small change in the T- and //-dependence 
of hnn, the coefficient that does not depend at all on Qp. Nevertheless in both analyses, fe^n 
and hncf) have the same behavior of initially falling very sharply from positive to negative 
values as the chemical potential n is increased from = 0, and gradually becoming zero 
as — > 4t. h(^(^ falls much faster with x in the exact computation than in the previous 
one, but still more gradually than hnn or ^n^- Again, this can be understood analytically 
in the low-tcmpcrature limit using (28)-(29), which gives bnn = bn(f> = and 6<^<^ = 
(7C(3)/3^m/7r'^?7.^)A/4(///4t) for /i 3> The qualitative agreement between the two 

approaches gives us confidence that the shape of the phase diagram in Figure 2 is indeed 
a robust prediction of the model, in the one-loop approximation. 
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The exactly computed coefficients (a„— a^)/^ (soUd) and (a<^— a^)/K (dashed) as a function of m/^j 
for the same values of /3A as in Figures 1 and 3. 
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The exactly computed coefficients 6„„T^/k, 6„^T^/«; and h^^T^ /k as a function of m/A, for the same 
values of (3 A as in Figures 1, 3 and 6. 



5.3) SO{5)-breaking due to 



One interesting issue we can explore within the exact computation is the extent to 
which 5*0(5) symmetry is broken by the assumption of d-wave SC pairing versus s-wave 
AF pairing. In Fig. 6 one sees that, despite our tuning the normahzation of Qp so as to 
preserve SO (5) near the FS, |Aa„| > |Aa<^| even at zero doping. We could have tuned gp 
such as to force Aa^ = Aa<^ at some temperature Tq, but the equality would hold only for 
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T = Tq. In fact we find that the ratio Aa<^/Aa^ scales hnearly with ln(/?/4t): 



Attn 



= 0.20 + 0.064 ln(/3/4t). (50) 

/u=0 



For comparison, the breaking in the quartic couphngs can be similarly parameterized in 
terms of the ratios 

(-^J ^ 0.082 + 0.10 ln(/3/4t); (-^) ^ 0.27 + 0.084 ln(/3/4t). (51) 



5.4) Derivative Terms in the Free Energy 

Using the above techniques and the procedure discussed in §4.5, we can evaluate the 
lowest derivative terms in the quadratic part of Fi, 

Fi = c^(V0)2 + c„(Vn)2. (52) 

Defining k.' = /3'^t/ (Aivh)'^ , i — e/At, fi = ///4t, and going to the limit of large (3t, we can 
write the coefficients as 

" ' C * ^^Mm - - - ""'^*''') '''' 

= ^-^^ ((1 - - (1 - IAD' MW) ; (84) 

" C © " " " '*"'"'^'*') '''^ 

~0(l//?t)c^ (56) 

To obtain (54) we used the identity (28), which is however only valid for ii 3> 1//3. For 
ji < 1/(3, the logarithmic singularities in Ao and A/2 are integrated over to obtain result that 
is finite as — > 0. The identity (28) gives a vanishing result when applied to (55), meaning 
that Cn is of order 1/pt compared to c<^, and thus must be computed numerically to find 
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the leading behavior in the large pt limit. Again, these statements hold only for ^ 1//? 
due to the singular behavior of the integrand at £ = 0. For /j, < 1/(3, we expect that 
C(t> ~ Cni just as was found for the coefficients of the quartic terms. 

These derivative terms will not be of further use in the present preliminary study, 
but may be useful in computing the higher order quantum corrections which we leave for 
future investigation. 

6. Conclusions 

In summary, our calculations lead us to the following conclusions: 

• Microscopic SOih) Invariance: We find, by construction, that microscopic models can 
indeed give rise to 50(5) invariance which relates the d-wave superconducting and antifer- 
romagnetic phases. This symmetry can emerge accidentally at long wavelengths, and so be 
relevant for the effective potential which determines which phase is energetically preferred 
(in a homogeneous system). Being a long- wavelength effect, the symmetry can be easily 
hidden in treatments which include all scales together. 

• Features of the Phase Diagram: Using our model we have computed the Landau-Ginzburg 
free energy whose minimization determines the phase diagram in the temperature-doping 
plane. We find that the coefficients of the Landau-Ginzburg potential have the right qual- 
itative dependence on temperature and chemical potential to allow for antiferromagnetism 
at half filling, but which is destroyed as the system is doped. It can also describe a 
superconductor which is favoured only for a range of nonzero dopings. 

• Significance of the SO (5) Symmetry: The class of models we considered encompasses 
S0{5) symmetry as a special case among the most general possible coupling functions. We 
find that if the symmetry is broken in the microscopic theory, there is no special tendency 
for it to be restored or further broken in the effective long-wavelength theory. The phase 
diagram is only of the desired form if there is a certain level of S'0(5)-breaking, whose 
fundamental origin is unknown. In the specific model considered in §5, the symmetry 
appears to be badly broken, since fp and gp are quite different functions of momentum, 
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yet by tuning their relative sizes it appears possible to obtain the correct phase diagram. 
At present it is therefore not yet clear whether nature demands an approximate SO (5) 
symmetry as the correct way to understand the relation between the AF and SC phases in 
cuprates, or whether 5'0(5)-symmetric models simply fall into the same universality class 
as the correct theory. 
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